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Abstract 

We use computer simulations to test a simple idea for mapping between long-time self diffu- 
sivities obtained from molecular and Brownian dynamics. The strategy we explore is motivated 

o: 

by the behavior of fluids comprising particles that interact via inverse-power-law pair potentials, 
which serve as good reference models for dense atomic or colloidal materials. Based on our sim- 
^ | ulation data, we present an empirical expression that semi-quantitatively describes the "atomic" 
to "colloidal" diffusivity mapping for inverse-power-law fluids, but also for model complex fluids 

o : 

'— '■ with considerably softer (star-polymer, Gaussian-core, or Hertzian) interactions. As we show, 
the anomalous structural and dynamic properties of these latter ultrasoft systems pose problems 

00 ' 

[ f° r other strategies designed to relate Newtonian and Brownian dynamics of hard-sphere-like 



o 
o 



X 



particles. 
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Computer simulations and statistical mechanical theories have long served as invalu- 
able tools for understanding relaxation processes that occur in fluid systems that range 
from molecular liquids to complex suspensions of Brownian particles.™ Depending on the 
resolution of the interparticle interactions used to model these systems, some descriptions 
of the microscopic dynamics are more appropriate to adopt than others. Simplistically, 
Newtonian (i.e., classical molecular) dynamics (MD) is suitable when the particles of in- 
terest are described by force fields with molecular-scale resolution. Modeling the motion 
of suspended Brownian particles, on the other hand, typically calls for coarser interparticle 
forces and dynamics that approximately account for the effects of the solvent 

Although the short-time behavior of a given model depends sensitively on its micro- 
scopic dynamics,- there is evidence that under certain conditions-e.g., dense fluids near 
the glass transition-the qualitative long-time behavior is largely independent of those 
details.- 1 ^ With this in mind, it is perhaps not surprising that researchers often select the 
type of microscopic dynamics to use in simulations based on other considerations, such 
as computational efficiency. A common example is adopting classical MD to explore the 
long-time dynamic behavior of model complex fluids with coarse, effective interactions, ig- 
noring the kinetic role of the implicit "fast" degrees of freedom, e.g, solvent.— - — Although 
the qualitative trends provided by these MD simulations have been valuable, there is still 
the question of how to map between the long-time MD data of such simulations and that 
which would have been produced assuming a different type of microscopic dynamics. 

Here, we take a first step toward addressing this issue. Specifically, we present computer 
simulation results that test a simple heuristic approach for mapping between long-time self- 
diffusion coefficients obtained from Newtonian and Brownian (i.e., overdamped Langevin) 
dynamics (BD). The method follows from a seemingly naive hypothesis that "what mat- 
ters" in such a mapping can be deduced from the behavior of a fluid of particles that 
interact via an inverse-power-law (IPL) pair potential [V(x) = ex _At ], where x = r/a, r is 
the interparticle separation, e and a represent characteristic energy and length scales of 
the interaction, and /i determines the steepness of the IPL repulsion. The IPL model is 
a natural reference system for dense atomic or colloidal fluids, whose static and dynamic 
properties are dominated by the repulsive part of their interactions. Furthermore, one 
can show^ - — that a one-to-one relation must exist between the following dimensionless 
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representations of the long-time diffusivities of an IPL fluid associated with the two types 
of dynamics: -DmdP 1 ^ 3 a/ 771 / and Dbd/D . Here, p is the number density, m is the 
particle mass, k-g is the Boltzmann constant, and T is temperature. -Dmd (Dbd) represents 
the long-time self diffusivity obtained from MD (BD) trajectories, respectively, and D is 
the value of Dbd in the dilute (p — > 0) limit. Based on previous work,— — it is clear that 
the relationship between .DmdP 1 ^ 3 V 'm / 'k-^T and Dbd/D for IPL fluids is approximately 
independent of \i for /x > 4. Hence, a more precise statement of our aforementioned hy- 
pothesis is that the quasi-universal IPL mapping relationship can also be used to estimate 
.Dmd from D> B d (or vice versa) for other types of fluids, perhaps including those with very 
different types of interactions and physical properties. 

To test this hypothesis, we carry out MD and BD simulations that probe the long-time 
dynamics of fluids of particles that interact via IPL potentials with fi = 8, 10, 12, 18, 
and 36. We also explore the behavior of fluids with particles interacting via ultrasoft 
Gaussian-core 25 [V(x) = eexp{— x 2 }], Hertzian 1 ^ [V(x) = e(l — a;) 5 / 2 ], and effective star- 
polymer— [V(x) = A{— \nx + B} for x < 1, and V(x) = ABx" 1 exp{(i? _1 — 1)(1 — x)} 
for x > 1] potentials. Here, B = / 1/2 /2, A = (20/9)k B T(B- 1 - l) 3 , and / is the star 
polymer arm number. The latter three model systems have received considerable attention 
recently in the theoretical soft matter literature,-^" 2 ^ and represent stringent test cases 
for our approach (and others) due to their distinctive dynamical trends; e.g., each exhibits 
a wide range of conditions where self diffusivity anomalously increases with increasing 
particle density, as opposed to the behavior of IPL fluids. The MD and BD simulations 
that we present here cover much of the computationally accessible phase space of these 
model fluids, including both equilibrium and moderately supercooled conditions (~ 10 3 
state points in total). As we show below, a potentially useful outcome of our analysis of 
this extensive data set is an empirical analytical equation that semi-quantitatively relates 
-Dmd and D>bd for these systems. 

The MD simulations of our study generate dynamic trajectories by solving Newton's 
equation of motion using the velocity-Verlet algorithm in the microcanonical ensemble-^ 
For IPL, Hertzian, and star-polymer fluids, they contain iV particles (N = 1000, 4000, and 
4000, respectively) and use integration time steps of 0.00005 ^/ma 2 /e , O.OIa/twt 2 Je, and 
0.001 ^/ ma 2 / k&T, respectively. The BD simulations presented here generate trajectories 
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by solving the Langevin equation in the high-friction limit using the conventional Brow- 
nian (Ermak) algorithm.- 1 ^ For the star-polymer fluid, they contain N = 4000 particles, 
and use a time step of At = O.Itb, where tb = mD /k-QT and D = 0.001 ^Jo^hoTjm. All 
simulations use a periodically replicated cubic cell with reduced volume V j a 3 determined 
by the number density pa 3 = Na 3 /V of interest. For the IPL, Hertzian, and star-polymer 
systems, the pair potentials are truncated at r cut = 1.4 — 3.9a (depending on p)^, la, 
and 3.6a respectively. Long-time tracer diffusivities from both MD and BD simulation 
trajectories are computed from the average mean-squared particle displacements via the 
Einstein relation. For the analysis presented below, we also include some previously re- 
ported MD simulation data 15 for the Gaussian-core fluid and BD simulation data^ for 
the IPL, Gaussian- core, and Hertzian fluids. The methods used in those studies are the 
same as those described above, and the required simulation parameters can be found in 
the original papers. 

Before examining the simulation data of the various model fluids discussed above, we 
first consider what should generally be expected about the relationship between the dimen- 
sionless diffusivities, D B d/D and -Drmd = Dymp 1 / 3 ^ m/ k B T . For example, to leading 
order in p, we know that 1 — D-qd/Dq oc p and -D^ D oc p, which together imply that 
1 - D BD /D oc in this limit. Furthermore, there is also evidence^ suggesting that 

-Dbd oc Dud for supercooled liquids near the glass transition (Dmd, -Dbd — > 0). Since 
Dud necessarily shows pronounced variations with small changes in p or T under these 
latter conditions, we also have Dbd/D oc -Drmd- The following expression, 

1 - D BD /D = (1 + Ci D RM d + C2 J Dr / md )~ 1 , (1) 

is an example of a simple heuristic functional form that interpolates between the aforemen- 
tioned characteristic "fast" and "slow" limiting behaviors. Below, we examine how well 
eq. [1] can describe the simulation data for a variety of fluids comprising hard to ultrasoft 
particles if c\ and C2 are treated as constants. 

Computer simulation data of D RM d plotted versus 1 — D B d/D for the IPL, Gaussian- 
core, and Hertzian fluids are presented in Fig. [TJ The data span the k B T /e — pa 3 plane 
(details in the caption), characterizing the relationship between MD and BD long-time 
diffusivities for these fluids in their equilibrium and moderately supercooled states. Also 
presented in Fig. [T]is a least-squares fit of the data using the form provided by eq. [Q As can 
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FIG. 1: Reduced long-time diffusivity from MD simulations -DmdP 1 ^ 771 /^^) 1 ^ 2 versus the 
fractional reduction in the long-time diffusivity from BD simulations (relative to the dilute value) , 
1 — D-q^/Dq. The blue circles correspond to the Hertzian fluid at temperatures ranging from 
0.002 — 0.02 e//cB and densities that range from 0.2 — 8.0 <r -3 . The red circles correspond to 
the Gaussian-core fluid at temperatures ranging from 0.004 — 0.2 e/k-Q and densities that range 
from 0.01 - 1.0 cr~ 3 . The green circles correspond to IPL fluids with exponents, fi, ranging from 
8 to 36 and values of the parameter pa 3 (e/kBT) 3 ^ that span the respective equilibrium fluid 
phases. The curve is a least-squares fit to the data for the three systems, and is given by eq. [Q 
with ci = 3.3176 and c 2 = 2.6645. 

be seen-despite the distinctive non-monotonic dynamic trends of the Gaussian-core and 
Hertzian fluids as a function of density^ 1 ^^^*^'^ -the data qualitatively behave as the 
equation predicts. In fact, as is illustrated in Fig. Slf, more than 98% of the simulation 
data for -Dbd for each of these model fluids are within 20% of the eq. [^estimation based on 
the simulated -Dmd- Furthermore, to the extent that the data shown in Fig. [1] reflect quasi- 
universal behavior, it suggests that the well-known, empirical dynamic freezing criterion 
for colloidal fluids (-Dbd/A) ~ 0.1),— has an analog in atomic systems (-Drmd ~ -03). 

As an illustration of how eq. [1] might be further used, we consider the star-polymer 
fluid 26 mentioned above. The soft, logarithmic repulsive interactions of this model are 
known to produce highly non-monotonic dynamic trends; e.g., diffusivity show two minima 
as a function of pa 3 (see Fig. [2] of this paper and Fig. 1 of Foffi et al.— ). When a model 
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FIG. 2: Reduced long-time self diffusivity [5md(^ 2 ^bT '/to) -1 / 2 and D-qd/Dq] for the star- 
polymer fluid plotted versus reduced density pa 3 from MD and BD simulations, respectively 
Fluids of stars with different arm numbers ranging from / = 25 — 52 (top to bottom) and re- 
duced densities pa 3 = — 2.6 are shown, (a) Results from MD simulations, (b) Results from 
BD simulations (symbols) and estimates (curves) based on simply substituting the -Dmd data of 
panel (a) into eq. Q] with constants given in the caption of Fig [TJ 

fluid displays such nontrivial behavior with one type of microscopic dynamics (e.g., MD), 
it is not obvious a priori whether the trends will necessarily be reflected when another 
type (e.g., BD) is employed. In fact, an earlier investigation of this system— made a point 
to report dynamic results from both types of simulations. What Fig. |2]d illustrates that 
is that one can, to a very good approximation, predict -Dbd/A) for this system by simply 
substituting its -Dmd data of panel Fig. [2^l into eq. [I] with no adjustable parameters (i.e., 
using ci and C2 from data in Fig. (U which did not include the star-polymer fluid). As is 
illustrated in Fig. Slf, the predicted values for 88% of the state points are within 20% of 
the simulation results. 

We are not aware of another mapping approach that can make predictions of simi- 
lar accuracy for both hard and ultrasoft particle fluids. One alternative strategy,— r — 
hypothesizes that P-Dmd/[p-D]o = -Dbd/A), where [pD] = ]im. p -> {pDyiv} . Although 
this relationship approximately holds for simple fluids with steep repulsions (the so-called 
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FIG. 3: Fractional reduction in the long-time diffusivity from BD simulations (relative to the 
dilute value), 1 — Dbd/Dq [panels (a)-(d)], and reduced long-time diffusivity from MD simula- 
tions DuDP^im/kBT) 1 / 2 [panels (e)-(h)], plotted versus negative excess entropy per particle, 
— s cx //cb- Dynamic data corresponds to that of Fig. [T] and [21 and — s ex data were computed via 
the free-energy-based simulation techniques discussed in the text. 



hard-sphere dynamic universality class), 35 we show that it breaks down qualitatively for 
fluids with ultrasoft interactions. Specifically, Fig. S2I illustrates that predictions based 
on this hypothesis for the star-polymer system are generally very inaccurate-except for a 
narrow region of phase space with extremely low pa 3 and high /-where particle overlaps 
are avoided. 

Note that a quantitative link between -Drmd and D-qh/Do is easy to establish for models 
where both can be expressed as single- valued functions of the same static quantity. For ex- 
ample, one can show^i that, for an IPL fluid, -Drmd and Dbd/D q are strictly single- valued 
functions of excess entropy (relative to ideal gas), s ex . In fact, the same is approximately 
true for other simple liquids that are "strongly correlating" and mimic a variety of static 
and dynamic properties of IPL systems.— 

Do Gaussian-core, Hertzian, and star-polymer fluids show excess-entropy scaling be- 
haviors similar to the IPL fluids? To check this, we compute s ex for these models using 
free-energy-based simulation methods. Specifically, we determine the density dependence 
of the Helmholtz free energy at high temperature using grand canonical transition-matrix 
Monte Carlo simulation.— We then carry out canonical temperature-expanded ensemble 
simulations^ with a transition-matrix Monte Carlo algorithm 38 to calculate the change in 



7 



Helmholtz free energy with temperature at constant density. Together, these simulations 
provide the excess Helmholtz free energy and excess energy, and hence s ex . Additional 
details on these simulations can be found elsewhere in our earlier papers.— 1 ^ 

The excess entropy scaling behaviors of 1 — -Drd/ Do an d -Drmd are plotted in Fig. [3] for 
all model fluids and state points shown in Fig. [Hand ED The main point is that, in stark 
contrast to the behavior of the IPL fluids, 1 — D-qd/Dq and -Drmd of the ultrasoft fluids 
are not (even approximately) single- valued functions of s ex . Hence, the success of the 
IPL-motivated mapping strategy between MD and BD diffusivities reported here cannot 
be explained by appealing to arguments about strongly-correlating fluids. 

As a final note, we emphasize that this mapping has yet to be tested for systems with 
structural and dynamic properties that are strongly influenced by attractive interactions, 
a class of fluids that we plan to investigate in the near future. 
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Figure SI. Ratio of long-time BD diffusivity estimated from eq. 1 in the text, -DBD,Pred, 
to that obtained from simulation, -Dbd, plotted as a function of D BD /D . The dashed red 
lines represent a 20% deviation of the predicted diffusivity from the value measured in 
simulation. 
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Figure S2. Long-time diffusivity of the star-polymer system from MD and BD simula- 
tions. Symbols are the same as in Figure 2 of the main text. Curves in panel (b) show 
prediction based on D BD /D = (D MD ) / (\pD] Q /p), where [pD} = lim p ^ P-DmeA 
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